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Abstract 

While interference alignment schemes have been employed to realize the full multiplex- 
ing gain of _ftT-user interference channels, the analyses performed so far have predominantly 
focused on the case when global channel knowledge is available at each node of the net- 
work. This paper considers the problem where each receiver knows its channels from all 
the transmitters and feeds back this information using a limited number of bits to all other 
terminals. In particular, channel quantization over the composite Grassmann manifold is 
proposed and analyzed. It is shown, for K-user multiple-input, multiple-output (MIMO) 
interference channels, that when the transmitters use an interference alignment strategy as 
if the quantized channel estimates obtained via this limited feedback are perfect, the full 
sum degrees of freedom of the interference channel can be achieved as long as the feedback 
bit rate scales sufficiently fast with the signal-to-noise ratio. Moreover, this is only one ex- 
treme point of a continuous tradeoff between achievable degrees of freedom region and user 
feedback rate scalings which are allowed to be non-identical. It is seen that a slower scaling 
of feedback rate for any one user leads to commensurately fewer degrees of freedom for that 
user alone. 

Index Terms 

Composite Grassmann manifold, finite-rate feedback, interference alignment, interfer- 
ence channel, MIMO, quantization. 
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I. Introduction 

The importance of the role played by interference management in wireless networks has 
prompted many researchers to analyze the interference channel from an information-theoretic 
perspective. While the capacity region remains unknown, many insightful characterizations 
have been developed through outer bounds [1], [2] and approximations [3]. The conventional 
wisdom was that if the strength of the interference is comparable to the actual signal, orthog- 
onalizing the signaling dimensions between users was best providing degrees of freedom 
per user in a A'-user interference channel. Abandoning this 'cake-cutting' approach however, 
[4] demonstrated the achievability of the sum degrees of freedom for this channel of y in 
this interference regime through what has come to be known as interference alignment (IA). 
This surprising result has spurred further research in this area, including in other aspects 
of the scheme as in [5] in analyzing limited multipath and variation with number of users 
in the system as well as in its applications as in [6] in mitigating inter-cell interference in 
a cellular network. In the specific setup of frequency- selective channels, [7] and [8] have 
suggested efficient schemes improving upon the original IA scheme in terms of data rates 
and multiplexing gains achievable through coding over a finite number of channel realiza- 
tions. Finally, while the notion of degrees of freedom is an academic construct, experimental 
verification of the benefits accruing from the interference alignment scheme is available from 
several papers, such as [9], where a combination of interference alignment and cancellation 
is shown to provide a significant throughput increase. 

The requirement of perfect channel state information at the transmitters (CSIT) by the 
IA scheme in [4] is, of course, practically unrealizable for a time- variant or frequency-selective 
system and this issue has recently begun to receive considerable attention ([6], [10], [11]). 
In particular, [6] analyzed the impact of imperfect channel knowledge of the sum mutual 



information achieved by interference alignment, when applied to the downlink of a cellular 
network using Orthogonal Frequency Division Multiplexing Access (OFDMA). On the other 
hand, by noting that perfect knowledge of the channel realization is a reasonable assumption 
to make for a static channel, [10] analyzed a M x M Gaussian interference channel with real 
or complex static coefficients and employed a result from Diophantine approximation theory 
to prove the achievability of ^M- degrees of freedom. 

It has been shown recently - for a frequency- selective single-input single-output (SISO) 
setup - by [11] that the full spatial multiplexing gain of y can be obtained even under 
conditions of limited feedback as long as the feedback rate exceeds K[L — l)logP bits per 
receiver, where L is the number of taps in the channel between any pair of nodes and P is the 
total power available with the transmitting sources. Their scheme involved a Grassmannian 
line quantization - as in [12] - of the channel vectors by the channel-aware receivers and 
their feedback using a limited number of bits broadcast to all the nodes in the network. 
The transmit and receiver beamforming vectors are then calculated at each node by treating 
these channel estimates as the actual channel realizations in the original IA procedure. While 
such an implementation would not align the interference perfectly, [11] demonstrated that 
their scheme keeps the interference power bounded in the relevant-signal dimensions, thereby 
achieving the full multiplexing gain. Further, an analogous result is provided by the same 
authors in [13] for the MISO network with the number of antennas at each transmitter 
being more than or equal to the number of users. Note that in such a scenario, interference 
alignment is not needed, as zero forcing combined with beamforming suffices to attain the 
maximum sum degrees of freedom. 

The sum degrees of freedom for the 2-user multiple-input multiple-output (MIMO) in- 
terference channel - with (Mi,M 2 ) antennas at the two transmitters and (Ni,N 2 ) antennas 



at the two receivers - was carried out in [14], where the achievability depended solely on 
beamforming and zero-forcing techniques. An extension to the K-user MIMO case with M t 
antennas at each transmitter and M r antennas at each receiver was provided in [15]. The 
sum degrees of freedom was bounded within close lower and upper bounds which coincided 
whenever the ratio ^|^'^ r | was an integer. In particular, this means that the spatial 
multiplexing gain is precisely known for both the single-input multiple-output (SIMO) and 
the multiple-input single-output (MISO) cases. The lower bounds obtained in [15] depend 
on IA with perfect channel knowledge at each node. 

In this paper, we explore all these multiple antenna cases, for a frequency-selective setup, 
under the regime of limited feedback. In the SIMO case with R antennas at each receiver, 
we find that as long as each receiver transmits no less than Nf = K(RL — 1) logP bits using 
quantization based on the composite Grassmann manifold, the nodes can utilize these channel 
estimates within the ambit of the IA scheme of [15] and still achieve the complete spatial mul- 
tiplexing gain. The MISO case then follows immediately as a consequence of the reciprocity 
of alignment principle enunciated in [16]. Extending the same idea to the if-user MIMO case 
with M t antennas at each transmitter and M r antennas at each receiver, we find that our 
scheme of IA involving limited feedback attains the same degrees of freedom as the original 
IA scheme as long as each receiver uses no less than Nf = mm{M t ,M r } 2 K(RL — l)logP 
bits, where R = gfffiffi ] • In each codebook over the composite Grassmann mani- 

fold is employed to jointly quantize the normalized channel directions. We further generalize 
these results by proving that if each user transmitted a fraction of the afore-mentioned num- 
ber of feedback bits, then the user can still attain a proportionate fraction of the degrees 
of freedom promised above. In particular, if user i fed back ctj • Nf (0 < a-i < 1) bits, then 
user i achieves degrees of freedom equal to a times what that user would achieve with global 



perfect CSIT. Interestingly therefore, a slower scaling of feedback rate for any one user leads 
to commensurately fewer degrees of freedom for that user alone. 

This paper is organized into five sections. Section I introduces the topic and delineates 
the notational convention followed in this paper. Section II computes the feedback scaling 
rate in the case of SIMO and MISO interference channels. Section III explores the MIMO 
M t x M r case and shows that interference alignment with limited feedback can attain the 
best-known lower bound on the sum degrees of freedom for this network, which coincides with 
the actual sum degrees of freedom whenever max j- ^'^f/l is an integer. Section IV evaluates 

° mm{Mt,M r \ ° 

the impact of a reduced feedback rate on the achievable degrees of freedom of the system. 
This section also interprets our result in light of some earlier results in this direction and 
discusses insights drawn from our analysis. Section V concludes the paper. 

Notation: K. and C represent the real and complex fields, respectively. If z e C, then z c 
represents its complex conjugate. The superscripts H and * represent the hermitian conjugate 
and the transpose of a matrix, respectively. Extending the conjugate notation to vectors, we 
denote as v c = [v^v^, ■ ■ ■ ,v°Y when v — [v±, v 2 , ■ ■ ■ ,f n ]*- The square brackets [.] and the letter 
I denote the time-slot index, usually varying from to L — 1. The circular brackets (.) and 
the letter m denote the frequency-slot index, usually varying from to N — 1. The symbols 
o and <S> represent the Hadamard and Kronecker products of matrices, respectively. We use 
the term 'const' to denote any constant independent of the power P, whose value might 
change from equation to equation. This notational device greatly simplifies the appearance 
of many equations and enables us to concentrate on the relevant portions of the calculations. 
CN represents the circularly symmetric complex normal distribution. All logarithms in this 
paper are taken with respect to base 2. 



II. For Single-Input Multiple-Output (SIMO) Systems 
A. System Model 

There are K single-antenna sources Si, S2, ■ ■ ■ , Sk with a message each for their respective 
single-antenna destinations Di,D 2) ■ ■ -,D K . In our K-user SIMO interference channel model, 
we assume that each transmitter has one antenna and each receiver has R antennas. The 
L-tap response between the source Sk and the destination Di is given by h ijk [l] G C Rxl , I G 
{0,1,..., L — 1}. We shall assume that these coefficients are drawn i.i.d from a continuous 
distribution such that these values are bounded with probability one and that the channel 
remains in the outage setting where the values h^l] do not change during the transmission of 
the signal. Let us define the matrix T ijk G C LxR with the L rows /i* fe [0], /i* fc [l], ■ • • , h\ k[L — 1] 
and denote its R columns as q^l],^^], . . . ,Cj jfc [i2]. 

( h\M \ 

= (ci tk [l],Ci tk [2],...,Ci tk [R]j ■ 

\K k [L-l]J 

The use of an OFDM-type cyclic signal model transforms this into iV parallel frequency 
flat channels, hi tk (r) G C Rxl , r G {0,1,..., TV — 1}. The input-output relations are described 
by a set of R equations, 

Vi(r) = h^x^r) + ^ /^(r)x fc (r) + Zl (r). (1) 

Here yi(r),Zi(r) G C Rxl are the channel output and the i.i.d. noise on the r-th tone at 
the i-th receiver, Xi(r) G C is the channel input on the r-th tone at the k-th transmitter and 
hi,k{ r ) ^ C Rxl is the channel vector from S^ to Di on the same tone. By defining the column 
vectors 

y^^O),^!),...,^-!^^ 1 , 



Ti _A 
i,k — 



[x i (0),x i (l),...,x i (N-l)]eC 



Nxl 



[ Zi (0), z t (l), Zi (N - 1)] ~ CN(0, N I) G C 



NRxl 



and the block diagonal matrix 



H iJk 4 diag{^( ), ^(1), . . . , h° fc (iV - 1)} G C 



NRxN 



we get the following equation 



Hi = H i4 X { + ^ H i,kXk + Zi. 



(2) 



Let us define the matrix F iyk e C NxR with the N rows h\ tk (0), h\ k (l), . . . , h\ k (N - 1) 



The receivers are assumed to have perfect knowledge of their respective channels, i.e. 
each destination D t knows T ijk e C LxR \/k G {1,2, . . . K} perfectly. As in the SISO case ana- 
lyzed in [11], we assume there exist error-free dedicated broadcast links from the destinations 
to every other node in the network. During an initial channel feedback phase, the receiver 
broadcasts its channel state information using Nf bits of feedback. This is followed by the 
data transmission phase. The maximal rate of communication between Si and Di such that 
the probability of error is driven to zero as the block-length goes to infinity is denoted as Ri 
with R sum = J2f =1 Ri- The sum degrees of freedom is defined as d sum = limp^oo ^fp-, where 
P is the total power constraint on the transmitting nodes. 




\hl k (N-l)J 



B. Hamming Bound on the Composite Grassmann Manifold 

The complex Grassmann manifold G n ^ viewed as a set, is the collection of all the k- 
dimensional subspaces within a n-dimensional Euclidean space C n . An alternate viewpoint, 
which is often useful, is to view it as an equivalence class of n x n unitary matrices defined 
in G ri: k — U(n)/(U(n — k)xU(k)). A point on G n ^ is then given by 



where Q G U(n). The study of the Grassmann manifold for wireless communication scenarios 
was boosted by the relation found between the manifold and non-coherent communication as 
reported in [17]. [12] utilized Grassmannian vector quantization in their analysis of finite-rate 
feedback MIMO point-to-point links. This line of thought has been explored in great detail 
by many authors, as in [18] and [19]. On the theoretical front, [20] computed the volume 
of a geodesic ball in G ny k and used the same to analyze MIMO systems under beamforming 
direction feedback. 

Multidimensional Grassmann analysis (i.e. on G n ^ when k ^ 1) cannot be directed in 
the current interference alignment scheme for a 3-user system, since if all that the users 
know is the subspace spanned by { p^jj ■> prfil ' ||fei3|| ^ ne y wom d not be able to align their 
vectors so that they are separable at receiver one. A more precise characterization in terms of 
the actual channel directions is indispensable here for the accomplishment of the full spatial 
multiplexing gain. Such a representation is precisely what is provided by procedures on the 
composite Grassmann manifold. 

In [21] the composite Grassmann manifold G™ k is formed by taking the direct sum of m 



[Q] 




Q k o 

Qn-k 




copies of the Grassmann manifold G n ^, i.e. 

m copies 

On the Grassmann manifold, a commonly used distance metric is the chordal distance d c . 
For the particular case of k — 1, it reduces to d 2 .(v 1: v 2 ) = 1 — \v^v 2 \ 2 - One can extend this 
distance to G™ x as follows: If P,Q G G™ 1: then P = [pi, . . . , p m ], and Q = [q 1: q m ], 
where pi, qi G G nji V % e {1,2, . . . ,m}. 

m 

d 2 (P,Q) 4 d 2 M, qi ). 
i=i 

In our analysis, we would need only G^ L1 ] so we particularize the succeeding discussion in 
this subsection to this particular manifold. 

One can define a ball of radius 5 around a point P G Gr L1 as 

B P (S) = { Q G Gr L1 | d(P,Q)<5}. 

Let Vol(A) denote the measure of the set A, with the measure being derived from the above 
defined distance metric. Since G^ L 1 is an homogenous space, Vo\(B P (S)) = Vo\(Bq(S)) V P, Q G 
Gft L1 . Hence, we can denote a ball of radius 5 in the composite Grassmann manifold as B(5) 
without any reference to its central point. Based on our distance metric, let us construct 
a maximal packing of spheres on the composite Grassmann manifold G^ L1 such that the 
minimum distance between the centers of any two spheres is more than S. This is the precise 
analogue of the Grassmannian sphere-packing problem considered in [22]. Let us choose this 
maximal packing code as our quantization codebook C. The number of codewords possible 
in such a maximal packing is given by the well-known Hamming bound, as follows: 

\C\ < 1 



Here, \C\ represents the number of codewords in our code C, and fj,(B(5)) represents 
the normalized volume of a geodesic ball in the composite Grassmann manifold, defined as 
fi(B(5)) = -^(g^) • ^ e no ^ e ^ na ^ smce the distances on the manifold are measured along 
geodesies (formed according to the choice of metric), the ball described below is formally 
called a geodesic ball in the manifold. 

A general result in [23] states that for any valid distance metric, the normalized volume 
of the ball is given by the following expansion : 

fi(B(S)) = const. 5 dimG «V(l+o(5 2 )). (3) 

Here, dimG^ L x is the real dimension of the composite Grassmann manifold and is given by 
KdimG RL:1 = 2K(RL-1). 

Substituting 2 N i = \C\, and noting that in the regime of increasing high Nf, the o(S 2 ) 
term can be ignored, we write 



N const const 



22K(RL-l) 

A realization x G Gr L1 shall be encoded using Nf bits corresponding to the index of the 
codeword in C closest to it, i.e. the quantized version of x shall be 

q(x) = argmin [d(x, Xi)]. 

Xiec 

For our code, produced from a maximal packing on the composite Grassmann manifold, the 
maximum distortion experienced by a realization is bounded by the minimum distance of 
the code as follows: 

. const 

A max < 5 < Nf , (4) 

22K(RL-1) 



where A max = max [d(x,q(x))}. 



Setting — to be equal to we get 

2 2K(RL-1) 



N f = K (RL — 1) logP bits. 



(5) 



C. Proposed Scheme 

The destination node Di, knowing the matrix 7$ & perfectly, forms a i?L-length norm- 



one vector v. 



i,k 



A vcc (Tj,k) 
~ ||vec(T iifc )|| 



. This calculation is performed for all k G {1,2, ... ,K}, creating 



a point on G^ L1 denoted by Qi = [v it i, v ij2 , ■ ■ ■ ,v^k\- This is quantized over our 2^-level 
codebook C on the composite Grassmann G^ L1 and is reconstructed by other nodes as 
Qi = [vi,i,Vi t 2, ■ ■ ■ ,Vi,k\- The i?L-length vector v ijk is re-organized - similar to T i)k - in form 
of a L x R matrix, 



Qi,k 



( KM \ 



(d it k[l],d itk [2], . . . ,d ijk [R]) ■ 



\K, k [L-l\) 

The nodes zero-pad the L-length vectors d i:k [m], me {1, 2, . . . , R} to iV-length and take 
their DFTs to obtain di }k [m}. These are arranged as shown to form a N x R matrix, 



Qi,k 



( KM \ 

KM 



= ( d i)Jfc [l],d i)fc [2],...,d. 



i,k I 



\wl k (N-l)J 
The nodes also form the RN x N matrix Wi k as 



W i%k = diag{^, fc (0), w itk (l), w itk (N - 1)} e C 



RNxN 



Now, we follow the IA scheme as in [15] by treating the above W ijk as the actual H i>k 



and utilize their interference alignment scheme of finding the direction vectors. Let 116N, 

T = KR{K-R- 1), and (6) 

N±(R + l)(n+l) T . (7) 
We shall be coding over this iV-symbols to achieve the following degrees of freedom : 




R{n+l) r ie{l,2,...,R+l}; 
R(n) r ie{R + 2,...,K}. 



Now the transmitter and receiver pair i choose di different transmit and receive beamforming 
vectors, viz. e C RNxl and t>™ G C Nxl , respectively. As long as L scales sufficient fast 
with T (cf. bandwidth scaling in [5]), one expects to find the vectors u™ an d v™ satisfying 
the following properties. 

\(uT) H W hZ (vT)\>c>0, 
Vie{l,2,...,K},me{l,2,...,d i }, (9) 

Vi E {1,2,..., K} and Vm^pG {1,2,...,^}, (10) 

«) H W hk (v P k ) = 0, 

\/i^ke{l,2,...,K} 
and Vme{l,2,...,d i },pe {1,2,. ..,4}- (n) 



The source Sk formulates dk independent symbols - x\, . . . ,x d k k G C, which are sent along 
the directions v\,... ,vf. k G C 7Vxl as follows : 



a* = X>* w ™' (12) 

m=l 



where \ \v™\ \ = 1 and -E(|x™| 2 ) = jfj-- Note that in the scheme of [15], 



i) m — ■ ■ ■ — v m 
v i — — v R+li 



&ndv p R+2 = ••• = < 



for all valid values of m and p. We have implicitly assumed here that the number of users 
K is greater than the number of receive antennas R per node in the SIMO system. In the 
case of K < R, no interference alignment is necessary as mere zero-forcing would attain the 
maximal attainable K degrees of freedom. 

D. Achievability Result 

Theorem II. 1: The interference alignment scheme delineated above for a general i^-user 
single-input multiple-output (SIMO) interference channel with one antenna at each trans- 
mitter and R antennas at each receiver, where each pair of nodes in the network has a L-tap 
frequency selective channel between them, achieves the full spatial multiplexing gain of 
as long as each destination transmits more than K(RL — l)logP bits of feedback, where P 
represents the total power available with the transmitting nodes of the network. 

Proof: We recall that using the Hamming bound on the composite Grassmann manifold 

G RL,1> 0Ile § etS 

const 

A max < Nf . (13) 

22K(RL-1) 

Setting — as -p, one obtained that Nf = K(RL — l)logP. Our aim is to show that 

2 2K(RL-1) 

this scaling of bits is sufficient to bound the interference terms in the rate expression. The 



destination Di projects the received signal y i onto the di directions given by u™ G C RNxl ,m G 
{1,2,. ..X}, 

+ ( 14 ) 

k^i p=l 

Let us define 1 R = [1, 1, . . . , 1]* G R Rxl , 

hk = [hUO),hl k (l),---,K k (N-l)}eC RN >< 1 , 

and w hk 4 K ifc (0)>* fc (l), . . . - 1)] G C™* 1 . 

Let us define (u™) c o (v p k ® l R ) as 6™/ G C RiVxl . Then, we can denote (u™) H H i:k v p k as 
h i k b^ k . Using this notation, 



dk 

fc^i p=l 



We again choose the input symbols x™ to be i.i.d Gaussian and since the receiver knows 
6™' m and h iti , it can treat the other interference as noise to obtain a rate of 

1 * / -^-\h H -b m ' m \ 2 \ 

A 2 = EE^rlO^I 2 - (is) 



is the interference (treated as noise) caused by transmitter % involving messages other 
than the one being currently decoded by receiver i. I i 2 is the interference caused by the 
transmitters other than transmitter i. 

Our three conditions on the vectors uT and vY 1 can also be re-written now as 



~Hi,m,m\ n \ j ■ 

WiAl I > c> V?, m, 



(19) 



We observe that hi^w^i and 6™' p are all i?iV-length vectors; and that 



(20) 
(21) 



N-l 



R 



Wi 



rn=l 



£||tM")ll 2 = £ll4,i["# 

n=0 



/"III 2 = 11-Oi.ill 2 = 1. 



m=l 



Since, one can always append vectors to the orthogonal vectors w iti and b^f to form a basis 



for C RN , it follows that, 



|^,i|| 2 > \{Ki) H w i:j !\ + 



(M 



im,p 
H i,i 



m,p I 
i,i I 



We thus get, 



< 



K<1, 1 

p - ndi^pvir-KM^r) 



Kd, 



Now, note that 



(hi,i) H _ ^(Ki) H j r , 

— = Wa = > — = 

IIM ' ^ IIM ' L J 

m =i IIMI 



The last line follows from the observation that, 



This leads to 



7V-1 



II 



m=l 



R 



£iim«)ii 2 =£ii* 

ra=0 m=l 

£ II c mNH 2 = ||vec(T i>fc 



m || 2 
i.i II 



p 


Km 2 






P 

" Kd % 


■ \Kf\ 


him 2 


(- 


(hi,i) H 

IIM 


P 

Kdi 


■ Wf\ 


i 2 iim 2 - 


(i- 




P 

Kdi 


■ m p \ 


l 2 IIM 2 - 


d 2 c {v iA , v iti ) 


P 

~ Kdi 


■ \Kf\ 


HIM 2 - 


K 






fc=i 




P 

Kdi 


■ m p \ 


i 2 iim 2 - 




ii Qi) 



J l,l 



< P. const. A, 



< P. 



const 



2k(rl-i) 

P. const. 



const. 



Hence, i^i is bounded. A similar argument using -jL=\h ik b™'£'\ shows that I^i is also 
bounded independent of power P. Further, as P — > oo (and the number of feedback bits 
Nf — > oo), the quantization error becomes vanishingly small by equation (4). This implies 
that, as in [11], that 



Now, the equation (16) yields the following result: 

r i . Rsum 

d sum = hm - 

P->00 lOg -T 

x * io g (i+ fi ; 1+f . 2+Aro 

i=l r, 



1 m=l 6 



N 

(R + l)fl(w + l) r + (K - R - l)Rn v 
~ (i?+l)(n + l) r 

The last step above follows from equations (6), (7) and (8). Taking the supremum over all 
values of the auxiliary parameter n, we obtain the desired result that d sum = which is 
the same as was obtained in the case of perfect channel state information at all the nodes in 
[15]. 

■ 

The MISO case follows immediately by application of the reciprocity of alignment scheme 
discussed in [16]. 

III. For Multiple- Input Multiple- Output (MIMO) Systems 

The precise degrees of freedom for the general M t x M r K-user interference channel 
have not yet been quantified. Gou and Jafar [15] established lower and upper bounds which 



coincide when R = I max f^ , fM I is an integer. We denote these as dH^ and dfjL below. 

L mm{iwt,Mr| J sum sum 



d 



UB 



K.min{M t ,M r } K < R; 
< - ) wyt y A' ... R . 



"'sum 



In the case of K < R, the <i^^ is obtained by noting that in the event of facing no interference, 

the d sum is upper bounded by K times the degrees of freedom of the single-user MIMO 

channel as calculated in [24]. In the case of K > R, a tighter upper bound is obtained by 

first reducing the K-user channel to different R + 1-user configurations; then by making the 

first R of the transmitters and the last R of the receivers cooperate; and finally by invoking 

the results presented in [14] on the 2-user case. 

The lower bounds are as follows: 

K.mm{M t ,M r } K < R; 
^.mm{M t ,M r }X K > R. 

The case of K < R is handled using just beamforming and zero forcing. The case of K > R 

is proved using interference alignment with perfect channel knowledge. We show that an 

appropriate scaling of the feedback bits suffices to achieve this spatial multiplexing gain. 

Theorem III.l: In a general i^-user interference channel where each transmitter has M t 
antennas, each receiver has M r antennas, K > R= [ ^^{mI'm^ ' anc ^ eacn P a ^ r °^ n °des in the 
network has a L-tap frequency selective channel between them, an interference alignment 
scheme under limited feedback achieves the same degrees of freedom as the interference 
alignment scheme with perfect channel state information as long as the receiver employs more 
than min{M t ,M r } 2 K(RL — l)logP bits of feedback, where P is the total power available 
with the transmitting nodes of the network. 

Proof: Due to the reciprocity of alignment idea [16], we can concentrate our attention 
on the case of M t < M r without any loss of generality. From each receiver node, we discard 



M r — RM t antennas. Now, by restricting the cooperation allowed amongst the transmitters 
and the receivers, we can treat this K-v&ei M t x RM t channel as a SIMO KM t -user 1 x 
R channel. As long as each user - utilizing composite Grassmann quantization - sends 
(KM t ).(RL — l).logP bits of feedback, we know from the preceding section that a spatial 
muliplexing gain of can be achieved. This value of matches with the desired inner 
bound achieved with perfect channel knowledge in [15]. Combining back the receivers to get 
the if-user channel, we conclude that each node needs to transmit at least M^K{RL — 1) 
bits of feedback. ■ 

IV. Remarks and Discussion 
A. Impact of smaller feedback scaling rate 

Through theorems II. 1 and III. 1 , we have established the achievement of the maximal 
sum degrees of freedom (d sum = ^ ' ), when the feedback by each user scales as some Nf 
bits. However, in an interference network, it may not be possible on part of each user to scale 
its feedback bits at this uniform rate. Obversely, it would be desirable to have a tradeoff 
between the provision of a flexible feedback rate for an individual user and the degrees of 
freedom obtained thereby, with the system-level strategy being optimized (as in Sections II 
and III) for achieving the Pareto-optimal point maximizing the sum degrees of freedom of 
the network. This tradeoff provides additional flexibility to the system designer in terms of 
allocation of resources for feedback and obtainable rates is given by the theorem below. 

We denote the degrees of freedom achieved by user % as di = ^, where di = limp^oo 
and Ri is the rate achieved by user % over N blocks. 

Theorem IV. 1: If the feedback rate employed by receiver i in the interference channel 
(both the SIMO/MISO model of Section II and the MIMO model of Section III) scales as 



(Xi.Nf, for some < a« < 1, then user i can achieve ctj.dj degrees of freedom by using 
interference alignment. 

Proof: We analyze below for a fT-user SIMO 1 x R channel; the corresponding MIMO 
results follow in a similar manner. Intuitively, the quantization error is proportional to 
P 1_Qi and this acts as a principal component of the interference faced by user i, leading to 
a reduction of (1 — ai)di degrees of freedom from user i. 

As in equation (22), the rate achieved by user i can be lower bounded as 

di 



m=l 

The interference terms can be expressed as 



m=l v * 7 
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dk 



P =i - J * 



If we denote b 2 max = max iifcitniP \\b™f\\ 2 , and h 2 max = max^ \\h i)k \\ 2 , then 



K 



< b m ax- h max . const. P. ^// ( 2 (r, (',.,,; 



k=l 

2/ 



= P. const. cr(Qi, Qi) 
< P. const. A* .■ . 



Here, A maxi is the maximum quantization error possible while using the codebook of 
receiver % with 2 aiN f representations on G^ L1 . As in equation (4), we can claim that 

const 

A maXj j ^ a~Nf 

22K(KL-1) 

. 9 const const 
A . < = 

^ max, l — a i K{RL-l)logP p a 
2 K(RL-l) 

This leads to, 



h,i + h,2 < const. P 
=>MG(Ri) > 1(^.(1 -a i )d i ) = ^ = a i 4 

where MG{x) denotes the multiplexing gain of the quantity x defined as MG{x) = limp^oo ^p- 



In particular, if a\ = a 2 = ... olk — ol (say) , then d sum = cxYld=i% = a -JTp[ m 
a K-user SIMO 1 x R channel. A similar result would also hold for the MIMO M t x M r 
channel. Also, note that while the rate achieved by user i would depend upon {(Xj}f =1 , the 
degrees of freedom depends only on ctj as long as the other a/s are positive. 

For the case of M-user Rayleigh faded broadcast channel with M antennas at the trans- 
mitter, [25] had analyzed the specific strategy of zero-forcing (ZF) precoding with random 
vector quantization or RVQ-generated quantization codebooks to show that a feedback of 
a times the optimal feedback rate leads to a times the maximum sum degrees of freedom. 
The above theorem presents a counterpart to that result for the ii'-user MIMO interference 
channel with several differences. We do not restrict the number of users to be equal to the 
number of antennas. The fading distribution is more general. No constraint is placed on 
the quantization codebook. Even when formulated using pseudo-beamforming vectors b™*, 
interference alignment differs widely from the simple ZF precoding schemes with the bj 1 ^ 's 
arising from involved functions of the channel matrices. 



B. Shrinkage of 'radius of uncertainty' 

The work in [26] analyzed the operation of a multi-antenna broadcast channel with the 
transmitter basing its beamforming vectors on imperfect channel estimates received from the 
receivers. Particularizing to the two-user case, they considered a scenario where the norm- 
one channel direction to user one is fixed at hi and the direction vector to user two can be 
either h 2a or h 2b . Irrespective of the transmission strategy adopted at the transmitter, [26] 
showed that the square of the chordal distance between h 2a and h 2b must shrink as 0(P _1 ), 
i.e 1 — \h 2a h 2b \ 2 = 0(P~ 1 ), for the system to achieve the same degrees of freedom as the 
perfect CSIT case. The achievability for the same was shown via a simple beamforming 
scheme. 

An analogous formulation of our problem would be to note that the rate achieved by 
user i, as expressed in equation (16), is 

Ri>I(Xi;Yi | {4}f) 

1 di I P \T H h m ' m \ 2 \ 

>lyiog(i + ^-\hffiTr\ 2 ) 

~ N ^ V Kali l ' 1 hl J 

m=l v 7 

- ivE lo g( J ^ + J ^ + iVo). (22) 

m=l 

Defining the multiplexing gain of a quantity x as limp^oo j^p, we note that the multiplexing 
gain of the first term equals the same degrees of freedom as available for the user in this 
channel in the perfect CSIT case. Further, we can show that the multiplexing gain of the 
second term is zero if both 1^ and l i 2 are shown to be constant independent of P. This 
condition is satisfied, as seen from our calculations in Section II. D, if A^ ax < p, which leads 
to Nf = K(RL — l)\ogP bits in equation (5). From Section II. C, we note that A^ ax = 



maXg i£C d 2 (Qi, Qi) , where Qi e G^ L1 and C is the code of cardinality 2 N f used to quantize 
Qi. Equivalently, one can re- write the necessary condition as 

W,,*) <^VAe{U..,4 (23) 

Recall that v^k is the norm-one estimate of the channel vector; and is its channel 
estimate provided by the receiver. As in the sufficiency part of the proof in [26], the v^k were 
treated as the actual channel directions in the formulation of the beamforming vectors (u™ 
and i>£) by the nodes. In our analysis, the role corresponding to the beamforming vectors of 
the broadcast channel are performed by the pseudo-beamforming vectors b™* formed by the 
Hadamard product of the transmit (t>£) and receive (u™) beamforming vectors. Our results 
interpreted in this manner act as a counterpart for the case of interference channels to the 
sufficiency results in [26] obtained for the broadcast channel. The complicated dependence of 
the pseudo-beamforming vectors b^ p on the fedback channel information precludes finding 
a direct counterpart to the necessity part of their proof. 

Note that in our case, v it k is uniformly distributed over a space G\rl,i with uncountably 
infinite points, i.e. it lies in a ball of uncertainty around the estimate v^. Equation (23) 
then shows that the shrinkage of its radius of uncertainty as 0(P _1 ) suffices to attain the 
same degrees of freedom as the perfect CSIT case. 

C. Remarks 

• The connection of the pre-log factor in the feedback scaling answer to the dimension of 
the manifold over which the quantization happens is intriguing. In schemes involving 
beamforming feedback of a single vector of M-length, the dimension of this norm-one 
vector is 2(M — 1). Correspondingly, we find that (M — 1) log P bits of feedback suffices 



to obtain ideal-like performance. In the interference alignment scheme explored above 
in the SIMO and MISO cases, we obtain exactly | dirnGf^ : as the pre-log factor. Note 
that the real dimension of the composite Grassmann manifold G^ L1 is 2K(RL — 1). 
This line of thought is further explored in Appendix A. 

• Since our analysis in Section II. D only required the first-order term in the volume 
expansion of a ball in G^, the expansion obtained in equation (3) was sufficient. 
The structure of G^, however, permits the evaluation of the precise volume of the 
ball as follows. 

Theorem IV. 2: The normalized volume of a ball of radius 5 in the composite Grass- 
mann manifold G^i is given by 

The proof is provided in Appendix B 

• It would be interesting to analyze interference alignment schemes under finite-rate 
feedback for time-selective channels. This paper's approach as also the analysis earlier 
by [11] cannot be directly applied to time- selective channels without violating the 
non-causality requirement on channel knowledge at the receivers. 

V. Conclusion 

Multi-user interference channels with single or multiple antennas at each node are ana- 
lyzed in a frequency-selective setup, wherein the receivers with knowledge of their respective 
channels quantize the channel directions using a code-book on the Composite Grassmann 
manifold and broadcast them to all other nodes at a rate that scales sufficiently fast with 
the power constraint on the nodes. It is shown that an interference alignment scheme based 
on treating these channel estimates as being perfect is sufficient to attain the same sum 



degrees of freedom as the interference alignment implementation utilizing perfect channel 
state information at all the nodes. We also demonstrate a continuous tradeoff whereby an 
individual user can opt for a slower scaling of feedback bits and obtain proportionally lower 
degrees of freedom. 

Appendices 

A. Pre-log Factor in Feedback Scaling Rate 

In this appendix, we give an intuitive explanation for why the pre-log factor in the 
feedback scaling rate required for emulating perfect CSI performance involves a term corre- 
sponding to the dimension of a manifold. The performance measure seems to be immaterial 
here, viz. it can be the maximal spatial multiplexing gain, the capacity of a channel or 
the probability of error incurred in a scheme. One first isolates that channel parameter the 
knowledge of which enables the scheme (using operations at both the transmitter and the 
receiver) to attain its 'ideal' value. This channel parameter can be, for example, 

• a single beamforming vector h G C Lxl , \\h\\ — 1 =>- dimension = L — 1, 

• an input covariance matrix Q e C nxn , Tr Q < p =>- dimension = n 2 , 

• a set of directions W G Gr L 1 =>- dimension = 2RK(L — 1). 

Let us denote this channel parameter as q, and denote its substratum on which it takes 
values as the manifold M. If the real dimension of M is given by dimM = N, then we know 
that around each point on the manifold, we can establish a local coordinate system with N 
variables. For analysis, consider such a system around the optimal value of q denoted as 
q opt . If the receiver only has knowledge of q, it can decide to feedback this information to the 
transmitter using a finite number, say Nf, of bits. By following the argument we made in 
Section II, one can bound the maximum distortion sufferable by a point under quantization 



using a 2 N f -level codebook on the manifold as: 

const 

A max ^ N . 

2n 

Here, the codebook is the one that solves the packing problem of spheres for a given code 
minimum distance for the given manifold M. The bound then follows using the Hamming 
upper bound and noting a general formula for ball volumes in a general manifold. We note 

that if q varies uniformly on the manifold M, then the distortion-rate function - which is 

N f 

the minimum distortion over the choice of all the codes - also varies as the inverse of 2t ) 
as seen in [27] and [20]. 

The receiver in attempting to convey q opt using Nf bits causes an error, and the trans- 
mitter is able to reconstruct the same as some q. In the regime of high Nf, we can take 
q within the same local coordinate system that covers q opt . Isolating the key part of the 
performance measure affected by the inaccuracy in q as f(q), we get an appropriate Taylor 
expansion: 

/(?) « f(q op t) + (V/)'(A ? ). 

V/ represents the gradient of the function / and Ag represents the displacement of q from 
q opt . Note that this is a coarse approximation chosen primarily for illustration; to make it 
rigorous, one would have to choose a normal coordinate system based at q opt and calculate 
distances along radial geodesies emanating from it. 

Our intention is often to bound the difference between f(q) and f(q op t) by a constant 
independent of P; as we expect to see both f(q) and f(q pt) increase monotonically with P. 
The variation of the (V/) with P determines the behavior required of the Aq term. If (V/) 



varies as VP, then we need Aq to vary as -^=. This is ensured by 

II Ag f < i A 2 nax < ^ 

const 1 N 

— 77 — < - ^> iVf si — log _r as -T increases. 

This provides an intuitive reason for our results on the feedback scaling rate in this paper 
as well as other similar results in [11] and [25]. 

B. Precise Ball Volume in the Composite Grassmann Manifold 

We are interested in finding the precise expression for the normalized volume of a ball in 
the Composite Grassmann manifold G% v For the general case of G njP , p G {1,2, . . . ,n}, [20] 
computed the first two terms in a series expansion for the normalized volume of a ball. The 
calculation has been extended in [21] to the case of the Composite Grassmann G% p . In the 
special case of G Ut i, the ball volume expression, somewhat surprisingly, reduces to a single 
term, as is extracted below from a simple manipulation of the results in [12]. The result of 
[20] for the case of general p does not reduce by setting p — 1 to the expression in [12] due to 
a mistake in their Corollary 1, wherein they claim the equivalence of the series expansion to 
its first term by showing merely that the second term in the expansion vanishes in certain 
cases. 

The paper [12] considers a spherical cap defined by 

S l ( 1 ) = {heC t \ |H| 2 = 7 , \<KC l >\ 2 > 7o }. 

Here Cj is an unit norm vector. The area of this spherical cap is then calculated as 

2tt*V7(7 - To)*' 1 



^(7)) = 



(*-!)! 



By setting 7 = 1, t = n, 5 = y/1 — 70, the spherical cap can be seen to be a ball in G nj i of 
radius 5, 

S l ( 1 ) = {heC n \ \\h\\ 2 = l, d 2 (/*,a)<5 2 }. 



The volume of the ball of radius 5 is then seen to reduce to 

27r n£2(n-l) 

Vo\(B(5))-- ~ 



(n-1)! ' 

Under the chordal distance metric, [21] provides the volume of the entire Grassmann manifold 
G n ,i as ^-"ly. ■ ^his gi ves the normalized volume of the ball in G n> i as 

M5(5)) = 5 2(n - 1) , 

as compared to the expression in [20] which only yields fi(B(5)) = 5 2 ( n-1 )(l + 0.5 2 + o(5 2 )). 

To extend the exact result to the Composite Grassmann case, we fix a point P e G^ 1: 
i.e. P = [Pi, P 2 , . . . , Pk], Pi £ G n ,iVi £ {1, 2, . . . , K}. We vary a random variable Q — 
[Qi,Q2, ■■■ ,Qk] uniformly on G* x , viz. Q { ~ Unif(G f „ j i) Vi G {1,2, By defining 

random variables Xj = d 2 (Pi,Qi), we note that these are independent and identically 
distributed with the cumulative distribution function 



F Xi (x) = Pr{X t <x} = n{B{yft)) = x n -\ 

Since the random variable X { measures the square of the cosine of the angle between two 
vectors, it is bounded between and 1. We also get the probability density function of Xi 

as 

f Xi (x) = (n-l).x n - 2 .l M , 

with l[ 0) i] being the standard indicator function for the interval [0,1]. 

Let us define U = J2?=i X% and we claim that the density of U is given by 

Mx) = <„-!) v^J^L lm . 

This can be proved through induction. Assume that statement is true for Yld=i Xi, the 



density of U can be found through convolution of the densities of X K and Ylf=i -^i as 

JU{ ' T((K-l)(n-l)) 
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The corresponding cumulative density function Fu(x) is given for x G [0, 1] as 



Fu(x) = {n-l 



K T K {n-l) x K ^ 
' Y{K{n-l))K{n-l) 
T K (n) 



x K{n ^ 1) . 



T(K(n + 
If we define a ball in as 

B(8) = {QeG* 1 \d(Q,P)<S}, 
then its normalized volume is given by 

f i(B(5)) = Pr{U<5 2 } = F u (5 2 ) 

F K (n) £2K(n-l) 



T(K(n + 
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